
#ifndef _LES_SSF_H_
#define _LES_SSF_H_

#include <cmath>
#include "Tuna.hpp"

namespace Tuna {

namespace LES_SSF {


  /*!
   *  Turbulent viscosity is calculated here.
   */
    template<class T_sfield>
	void calcTurbulentViscosity(T_sfield &,
	                           const T_sfield &,
				   const T_sfield &,
				   const T_sfield &,
				   double, double, double);

  /*!
   * The process of selecting turbulent viscosity or not is done here.
   */
    template<class T_sfield>
	void selectTurbulentViscosity(T_sfield &,
                                      const T_sfield &, 
				      const T_sfield &, 
				      const T_sfield &,
				      double, double, double);
    

}

} // Tuna namespace

#include "../../src/LES/les_ssf.cpp"

#endif // _LES_SSF_H_ 
